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Abstract 



We evaluate the relic density of the lightest neutralino, the lightest supersymmetric 
particle, in the Minimal Supersymmetric extension of the Standard Model (MSSM). 
For the first time, we include all coannihilation processes between neutralinos and 
\^ , charginos for any neutralino mass and composition. We use the most sophisticated 

■ routines for integrating the cross sections and the Boltzmann equation. We properly 

' treat (sub) threshold and resonant annihilations. We also include one-loop corrections 

, to neutralino masses. We find that coannihilation processes are important not only 

' for light higgsino-like neutralinos, as pointed out before, but also for heavy higgsinos 

0^ ' and for mixed and gaugino-like neutralinos. Indeed, coannihilations should be included 

fH , whenever < 2|Mi|, independently of the neutralino composition. When ~ \Mi\, 

^ ' coannihilations can increase or decrease the relic density in and out of the cosmologically 

I interesting region. We find that there is still a window of light higgsino-like neutralinos 

Oh! that are viable dark matter candidates and that coannihilations shift the cosmological 

^ ' upper bound on the neutralino mass from 3 to 7 TeV. 



1 Introduction 

In the near future, it may become possible to constrain supersymmetry from high precision 
measurements of the cosmological parameters || , among which the dark matter density. 
It is therefore of great importance to calculate the relic density of the lightest neutralino as 
accurately as possible. 

The lightest neutralino is one of the most promising candidates for the dark matter in the 
Universe. It is believed to be the lightest stable supersymmetric particle in the Minimal 
Supersymmetric extension of the Standard Model (MSSM). It is a linear combination of the 
superpartners of the neutral gauge and Higgs bosons. 

The relic density of neutralinos in the MSSM has been calculated by several authors during 
the years @, ||, @, |[ § with various degrees of precision. A complete and precise calcu- 
lation including relativistic Boltzmann averaging, subthreshold and resonant annihilations, 
and coannihilation processes is the purpose of this paper. 

*E-maiI address: edsjo@teorfys.uu.se 
tE-mail address: gondolo@mppmu.mpg.de 
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As pointed out by Griest and Seckel Q one has to include coannihilations between the hghtest 
neutrahno and other supersymmetric particles heavier than the neutralino if they are close in 
mass. They considered coannihilations between the lightest neutralino and the squarks, which 
occur only accidentally when the squarks are only slightly heavier than the lightest neutralino. 
In contrast, Mizuta and Yamaguchi pointed out an unavoidable mass degeneracy that 
greatly affects the neutralino relic density: the degeneracy between the lightest and next- 
to-lightest neutralinos and the lightest chargino when the neutralino is higgsino-like. They 
considered coannihilations between the lightest neutralino and the lightest chargino, but 
only for neutralinos lighter then the W boson and only with an approximate relic density 
calculation. Moreover, they did not consider Higgs bosons in the final states. 

Drees and Nojiri |^ included coannihilations into their relic density calculation, but only 
between the lightest and next-to-lightest neutralinos. These coannihilations are not as im- 
portant as those studied by Mizuta and Yamaguchi. Recently, Drees et al. re-investigated 
the relic density of light higgsino-like neutralinos. They included coannihilations between the 
lightest and next-to-lightest neutralinos as well as those between the lightest neutralino and 
the lightest chargino. They do however only consider //, //' and jW^ final states through 
Z and W exchange respectively, and do not consider t- and u-channel annihilation or Higgs 
bosons in the final states. 

In this paper we perform a full calculation of the neutralino relic density for any neutralino 
mass and composition, including for the first time all coannihilations between neutralinos and 
charginos. We properly compute the thermal average, particularly in presence of thresholds 
and resonances in the annihilation cross sections. We include all two-body final states of 
neutralino-neutralino, neutralino-chargino and chargino-chargino annihilations. We leave 
coannihilations with squarks |^ for future work, since they only occur accidentally when the 
squarks happen to be close in mass to the lightest neutralino as opposed to the unavoidable 
mass degeneracy of the lightest two neutralinos and the lightest chargino for higgsino-like 
neutralinos. 

In section ^ we define the MSSM model we use and in section we describe how we generalize 
the Gondolo and Gelmini | pO[ formulas to solve the Boltzmann equation and perform the 
thermal averages when coannihilations are included. This is done in a very convenient way 
by introducing an effective invariant annihilation rate WcS- In section ^ we describe how we 
calculate all annihilation cross sections, and in section ^ we outline the numerical methods 
we use. We then discuss our survey of supersymmetric models in section ^ together with 
the experimental constraints we apply. We finally present our results on the neutralino relic 
density in section ^ and give some concluding remarks in section |[ 



2 Definition of the Supersymmetric model 

We work in the framework of the minimal supersymmetric extension of the standard model 
defined by, besides the particle content and gauge couplings required by supersymmetry, the 
superpotential (the notation used is similar to Ref. [Tl|| ) 

W = [~e*j,YEllHi ~ d^jYz^qliJ^ + u|,YyqlJ?^ - ^.HIhC) (1) 

and the soft supersymmetry-breaking potential 

^soft = {e*j,AEYEllHi + d^AcYcqli/^' - u^jAc^Y^qii?! - BfiHlH^ + h.c.) 
+Hl*mlHl + H'^*mlW2 
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+ ^MiBB + ^Ma (W^W^ + 2W+W-^ + ^M^gg. (2) 

Here i and j are SU(2) indices (ei2 = +1). The Yukawa couplings Y, the soft trihnear 
coupUngs A and the soft sfermion masses M are 3x3 matrices in generation space, e, 1, u, d 
and q are the superfields of the leptons and sleptons and of the quarks and squarks. A tilde 
indicates their respective scalar components. The L and R subscripts on the sfermion fields 
refer to the chirality of their fermionic superpartners. B, and are the fermionic 
superpartners of the SU(2) gauge fields and g is the gluino field, n is the higgsino mass 
parameter, Mi, M2 and M3 are the gaugino mass parameters, 5 is a soft bilinear coupling, 
and rrii 2 are Higgs mass parameters. 

For Ml and M2 we make the usual GUT assumptions 

Ml = tan^ 0w ~ O.SAfa (3) 

M2 = . 2^7 M3-O.3M3 (4) 

sm OwCis 

where a^-vi is the fine- structure constant and Us is the strong coupling constant. 

Electroweak symmetry breaking is caused by both H\ and acquiring vacuum expectation 
values, 

{Hl)=vi, {Hl)=V2, (5) 

with g^{v\ + V2) — 2m^, with the further assumption that vacuum expectation values 
of all other scalar fields (in particular, squarks and sleptons) vanish. This avoids color 
and/or charge breaking vacua. It is convenient to use expressions for the Z boson mass, 
"^z — 5(5^ + ff'^)!^! + ^2) ^'^d the ratio of vacuum expectation values tan/3 = vijvi. g and 
g' are the usual SU(2) and U(l) gauge coupling constants. 

When diagonalizing the mass matrix for the scalar Higgs fields, besides a charged and a 
neutral would-be Goldstone bosons which become the longitudinal polarizations of the 
and Z gauge bosons, one finds a neutral CP-odd Higgs boson A, two neutral CP-even Higgs 
bosons iJi,2 and a charged Higgs boson H^. Choosing as independent parameter the mass 
mA of the CP-odd Higgs boson, the masses of the other Higgs bosons are given by 

K42 _( rn\ cos^ 13 + rn\s\v? (3 + AM\i - sin/3cos/3(m^ -I- m|) -I- AA^^2 \ ia\ 



sin/3cos/3(m^ -I- rn|) -I- AA^gi ^^^111 /3 + ™z cos (3 + AA^ 



22 



(7) 



The quantities AAlf^- and A± are the leading log two- loop radiative corrections coming from 
virtual (s)top and (s)bottom loops, calculated within the effective potential approach given 
in (other references on radiative corrections are [Q). Diagonalization of Al|^ gives the 
two CP-even Higgs boson masses, niHi 21 and their mixing angle a (— 7r/2 < a < 0). 

The neutralinos are linear combinations of the superpartners of the neutral gauge bosons 
B, W3 and of the neutral higgsinos i?", i?" • this basis, their mass matrix is given by 



( Ml 


g'vi 





M2 

_gv2_ 
V2 



g "1 



_gv2_ 



(8) 
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where ^33 and S44 are the most important one-loop corrections. These can change the neu- 
trahno masses by a few GeV up or down and are only important when there is a severe mass 
degeneracy of the lightest neutralinos and/or charginos. The expressions for S33 and Sm are 

S33 = ~j^Y,^mbSm{29i)Re[BoiQ,bM)~Bo{Q,bM)] (9) 
S44 = -Y^Yt^mtsm{29i)Re[BoiQ,t,ii)^ BoiQ,t,i2)] (10) 

where nib and mt are the masses of the b and t quarks, lb — gmi, / V2mw cos (3 and Yt = 
gmt/V^mw sin P are the Yukawa couplings of the b and t quark, 6^ and 9j: are the mixing 
angles of the squark mass eigenstates {qi = q'LCOS^+ qnsmOq) and Bq is the two-point 
function for which we use the convention in Ref. ^, |lj] . Expressions for Bq can be found in 
e.g. Ref. [^. For the momentum scale Q we use as suggested in Ref. Note that the 
loop corrections depend on the mixing angles of the squarks which in turn depend on the 
soft supersymmetry breaking parameters Ajj and Ajj in Eq. (^) (or the parameters Ai, and 
At given below). 

The neutralino mass matrix, Eq. (||), can be diagonalized analytically to give four neutral 
Majorana states. 



X 



^ NaB + N,2W^ + iV,3^? + N,4Hl (11) 



the lightest of which, to be called Xi is then the candidate for the particle making up the 
dark matter in the Universe. The gaugino fraction, of neutralino i is then defined as 

^3 = l^'a|' + l^.2p (12) 

We will call the neutralino higgsino-like if Zg < 0.01, mixed if 0.01 < Zg < 0.99 and gaugino- 
like if Zg > 0.99, where Zg = Zg is the gaugino fraction of the lightest neutralino. Note that 
the boundaries for what we call gaugino-like and higgsino-like are somewhat arbitrary and 
may differ from other authors. 

The charginos are linear combinations of the charged gauge bosons and of the charged 
higgsinos , H2 ■ Their mass terms are given by 

{W- H^)M^± (^^+) +h.c. (13) 

Their mass matrix, 

M^. = (^'' '"0, (14) 
is diagonalized by the following linear combinations 

X" = UaW-+UaHr, (15) 
xt = VaW+ + V,2Ht- (16) 



We choose det(J7) = 1 and U*Mji±V'^ = diag(m~± , to-±) with non-negative chargino masses 
TOp^± > 0. We do not include any one-loop corrections to the chargino masses since they are 
negligible compared to the corrections ^33 and 844 introduced above for the neutralino masses 
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When discussing the squark mass matrix including mixing, it is convenient to choose a basis 
where the squarks are rotated in the same way as the corresponding quarks in the standard 
model. We follow the conventions of the particle data group [|l6j and put the mixing in the 
left-handed c?-quark fields, so that the definition of the Cabibbo-Kobayashi-Maskawa matrix 
is K = V1V2, where Vi (V2) rotates the interaction left-handed u-quark (d-quark) fields to 
mass eigenstates. For sleptons we choose an analogous basis, but due to the masslessness of 
neutrinos no analog of the CKM matrix appears. 

We then obtain the general 6 x 6 u- and d-squark mass matrices: 

= [^Q + mtm„ + Dl^l mt (A^, - /i* cot/3) \ 
" V (A[/ - Ai cot /3)m„ -f m„mt + D^j^l J ' ^ ' 

^.^/K^MlK + m^ml + Di^l m^(At, - ^* tan/3) \ 
V (Ai5-Mtan/3)m<j Ml + mlnid + D^j^l J ' 

and the general sneutrino and charged slepton masses 

Ml = Ml+Dl^l (19) 



2 ^ (mI + uiM + DIlI mt(Al,-/i*tan/3) 
(Ab -^tan/3)me M| + mtnie + ZJfj^jl 



(20) 



Here 

Dii^^m%co82P{T3f-efSm^ew), (21) 

= m| cos 2/3e / sin^ 0w (22) 

where T^f is the third component of the weak isospin and e/ is the charge in units of the 
absolute value of the electron charge, e. In the chosen basis, we have m„ = diag (tou, Wc m-t), 
rcid = diag(md,ms,mb) and nie = diag(me, m^, m^). 

The slepton and squark mass eigenstates fk {i>k with k = 1,2,3 and e^, Wfe and dk with 
= 1, . . . , 6) diagonalize the previous mass matrices and are related to the current sfermion 
eigenstates /lq and fjia (a = 1,2,3) via 
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ha = E/fc^FL, (23) 

fe=l 

fRa = Y.f'^^FR- (24) 

fe=l 

The squark and charged slepton mixing matrices Tul,r, Tdl,r. and Tel,r have dimension 
6x3, while the sneutrino mixing matrix Fj^i has dimension 3x3. 

For simplicity we make a simple Ansatz for the up-to-now arbitrary soft supersymmetry- 
breaking parameters: 

Act = diag(0,0,At) 

Ad = diag(0, 0, Ab) , ^ 

Ae - ^^^> 
Mq = Mu = Md = Me = Ml = mol. 
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This allows the squark mass matrices to be diagonalized analytically. For example, for the 
top squark one has, in terms of the top squark mixing angle di, 



UL 



p^ii 



smf 



(26) 



Notice that the Ansatz (25) implies the absence of tree-level flavor changing neutral currents 
in all sectors of the model. 



3 The Boltzmann equation and thermal averaging 

Griest and Seckel |^ have worked out the Boltzmann equation when coannihilations are 
included. We start by reviewing their expressions and then continue by rewriting them into 
a more convenient form that resembles the familiar case without coannihilations. This allows 
us to use similar expressions for calculating thermal averages and solving the Boltzmann 
equation whether coannihilations are included or not. 



3.1 Review of the Boltzmann equation with coannihilations 

Consider annihilation of N supersymmetric particles Xi (* = li • ■ • i with masses rrii and 
internal degrees of freedom (statistical weights) gi. Also assume that mi < TO2 < • • • < 
rriN-i < and that i?-parity is conserved. Note that for the mass of the lightest neutralino 
we will use the notation and mi interchangeably. 

The evolution of the number density Ui of particle i is 
dn ^ 

-E[r',(".-^r)-r,.(n,-n°^)]. (27) 

The first term on the right-hand side is the dilution due to the expansion of the Universe. 
H is the Hubble parameter. The second term describes XiXj annihilations, whose total 
annihilation cross section is 

<^ij = ^^). (28) 

X 

The third term describes Xi ^ Xj conversions by scattering off the cosmic thermal back- 
ground, 



^Xi] 



Y^aix^X^XjY) (29) 



being the inclusive scattering cross section. The last term accounts for Xi decays, with 
inclusive decay rates 

r^, = Y.^{x^^x,x). (30) 

X 
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In the previous expressions, X and Y are (sets of) standard model particles involved in the 
interactions, Vij is the 'relative velocity' defined by 



V _ (31) 



EiEj 



with Pi and Ei being the four-momentum and energy of particle i, and finally n°'^ is the 
equilibrium number density of particle Xi, 



9t 



(27r)3 



d'p./,: (32) 



where Pi is the three-momentum of particle i, and f i is its equilibrium distribution function. 
In the Maxwell-Boltzmann approximation it is given by 

h = e~^^/^. (33) 

The thermal average (aijVij) is defined with equilibrium distributions and is given by 

^"-"-^ ^ Id^P.d^p,fJ, 

Normally, the decay rate of supersymmetric particles Xi other than the lightest which is stable 
is much faster than the age of the universe. Since we have assumed i?-parity conservation, 
all of these particles decay into the lightest one. So its final abundance is simply described 
by the sum of the density of all supersymmetric particles. 



N 

n 



= Y.^i- (35) 



For n we get the following evolution equation 
dn ^ 

— = -SiJn - ^ (cTy Wy ) {n,nj - n^n^) (36) 

where the terms on the second and third lines in Eq. ( p7| ) cancel in the sum. 

The scattering rate of supersymmetric particles off particles in the thermal background is 
much faster than their annihilation rate, because the scattering cross sections cr'^^ - are of the 
same order of magnitude as the annihilation cross sections but the background particle 
density nx is much larger than each of the supersymmetric particle densities rii when the 
former are relativistic and the latter are non-relativistic, and so suppressed by a Boltzmann 
factor. In this case, the Xi distributions remain in thermal equilibrium, and in particular 
their ratios are equal to the equilibrium values, 

(37) 

We then get 
dn 



dt 



3Hn - {a,sv) (n^ - ) (38) 
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where 



oq „eq 
3 



= (39) 
3.2 Thermal averaging 

So far the reviewing. Now let's continue by reformulating the thermal averages into more 
convenient expressions. 



We rewrite Eq. (39) as 



(a^^v) = -2 = . (40) 



'^eq '^eq 



For the denominator we obtain, using Boltzmann statistics for fi, 

i i ^ i 

where K2 is the modified Bessel function of the second kind of order 2. 

The numerator is the total annihilation rate per unit volume at temperature T, 

ij ij 

It is convenient to cast it in a covariant form, 

Wij is the (unpolarized) annihilation rate per unit volume corresponding to the covariant 
normalization of 2E colliding particles per unit volume. Wij is a dimensionless Lorentz 
invariant, related to the (unpolarized) cross section through^ 

Wij = 'ipijVsaij = Aaij ^ (p^ ■ pjY - m^m^ = AEiEja^jVij. (44) 

Here 

[s - (m, + rrij)^] [s - (m, - rrij)^] ^'"^ 
P. (45) 

is the momentum of particle Xi (or Xj) the center-of-mass frame of the pair XiXj- 

Averaging over initial and summing over final internal states, the contribution to Wij of a 
general n-body final state is 

^^■^ internal d.o.t. f ' 



-fThe quantity Wij in Ref. W is Wij/i. 
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where Sf is a symmetry factor accounting for identical final state particles (if there are K sets 
of iVfe identical particles, k — 1, . . . , K , then Sj = Y[k=i -^fcO- particular, the contribution 
of a two-body final state can be written as 

""i^'-m^^i^ S /ixw^wr^n, (47) 

■' internal d.o.f. 

where pki is the final center-of-mass momentum, Sm is a symmetry factor equal to 2 for 
identical final particles and to 1 otherwise, and the integration is over the outgoing directions 
of one of the final particles. As usual, an average over initial internal degrees of freedom is 
performed. 

We now reduce the integral in the covariant expression for A, Eq. (|4^), from 6 dimensions to 
1. Using Boltzmann statistics for fi (a good approximation for T <m) 

A-V [ aaW- e-^^/^e-^^/'^—^ '^^^^ (i8) 

If' (2^)32i?, (2^)32i?, ' 

where Pi and pj are the three-momenta and Ei and Ej are the energies of the colliding 
particles. Following the procedure in Ref. pof we then rewrite the momentum volume element 
as 

cfpid^Pj = 4TT\p^\EidEi 4TT\pj\E-idEj ^dcosd (49) 

where 9 is the angle between and Pj. Then we change integration variables from E'i, Ej, 
6 to and s, given by 

E+ = Ei + Ej 

E^ = E,- Ej (50) 
+ 2EiEj 



+m'] + 2E,Ej - 2|p,| |pj | cos 61, 



whence the volume element becomes 

d^Pi d^pj _ 1 dE+dE-ds 
{2nf2E, {2TTf2Ej " (2^ 8 ' ^ ' 

and the integration region {Ei > mi, Ej > nij, \ cos9\ < 1} transforms into 

s > {ni, + mjf, (52) 
E^ > ^Ts, (53) 



E^ - E^ 



9 9 

m - — mj 



Notice now that the product of the equilibrium distribution functions depends only on E^ 
and not E^ due to the Maxwell-Boltzmann approximation, and that the invariant rate Wij 
depends only on s due to the neglect of final state statistical factors. Hence we can immedi- 
ately integrate over E^ , 



dE-=Ap.,j\j^ -. (55) 
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The volume element is now 



dE+ds (56) 



{2T:f2E, {2tt)32Ej {2tt)^ 2 
We now perform the E^ integration. We obtain 

^ = ^ E dsg,g,p,,W,,K, (^) (57) 



where Ki is the modified Bessel function of the second kind of order 1. 

We can take the sum inside the integral and define an effective annihilation rate Wefr through 

^gigjPijWij ^ gfpcsWcs (58) 

with 

1 



PoS = Pii = s - Amf. (59) 



In other words 



YPii 5? YV s{s-iml) g\ ' 

Because Wij^s) = for s < (jm + m^)^, the radicand is never negative. 
In terms of cross sections, this is equivalent to the definition 

sr^ Pij 9igj 



Eq. (57) then reads 

^ = ^ r dsp,sW,sKi ] (62) 



327r4 74„j2 \ T 

This can be written in a form more suitable for numerical integration by using pcff instead 
of s as integration variable. From Eq. (|59|), we have ds — Spcsdpcs, and 

^^^l dp.splsW.sKi (^^^ (63) 

with 



s 



Apis + (64) 



So we have succeeded in rewriting A as a 1-dimensional integral. 

From Eqs. (^3|) and (^), the thermal average of the effective cross section results 



/o°° dp.splsW.sKi (#) 



(65) 



10 
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m,o = 76.3 GeV 
m^o = 96.3 GeV 
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Figure 1: The effective invariant annihiliation rate Weff as a function of poff for model 1 
in Table ^. The final state threshold for annihilation into W^W~ and the coannihilation 
thresholds, as given by Eq. (|60|), are indicated. The X2X2 coannihilation threshold is too 
small to be seen. 



This expression is very similar to the case without coannihilations, the differences being the 
denominator and the replacement of the annihilation rate with the effective annihilation rate. 
In the absence of coannihilations, this expression correctly reduces to the formula in Gondolo 
and Gelmini pO[ |. 

The definition of an effective annihilation rate independent of temperature is a remarkable 
calculational advantage. As in the case without coannihilations, the effective annihilation 
rate can in fact be tabulated in advance, before taking the thermal average and solving the 
Boltzmann equation. 

In the effective annihilation rate, coannihilations appear as thresholds at ^/s equal to the 
sum of the masses of the coannihilating particles. We show an example in Fig. |l| where it is 
clearly seen that the coannihilation thresholds appear in the effective invariant rate just as 
final state thresholds do. For the same example, Fig. ^ shows the differential annihilation rate 
per unit volume dA/dpeS, the integrand in Eq. (|6^), as a function of PcS- We have chosen 
a temperature T = to^/20, a typical freeze-out temperature. The Boltzmann suppression 
contained in the exponential decay of Ki at high pcff is clearly visible. At higher temperatures 
the peak shifts to the right and at lower temperatures to the left. For the particular model 
shown in Figs. the relic density results fl^h^ — 0.030 when coannihilations are included 
and Op^/i^ = 0.18 when they are not. Coannihilations have lowered fl^h^ by a factor of 6. 

We end this section with a comment on the internal degrees of freedom gi. A neutralino 
is a Majorana fermion and has two internal degrees of freedom, g^o = 2. A chargino can 
be treated either as two separate species xt ^^'^ X7 : each with internal degrees of freedom 
9x+ ~ 9x~ ~ ^' simply, as a single species with g^± 

freedom. The effective annihilation rates involving charginos read 



4 internal degrees of 



^xfxf 



^?x- 



V i = 1, 



W + - 

x7x,- 



Xi X,- 



,4, j = l,2 



Vi,j = l,2 



(66) 
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10 20 30 40 50 60 70 80 90 100 



Peff [GeV] 

Figure 2: Total differential annihilation rate per unit volume dA/dpcs for the same model as 
in Fig. |l|, evaluated at a temperature T = to^/20, typical of freeze-out. Notice the Boltzmann 
suppression at high poff • 



3.3 Reformulation of the Boltzmann equation 



We now follow Gondolo and Gelmini |10 to put Eq. ( p38| ) in a more convenient form by 
considering the ratio of the number density to the entropy density, 

r = -. (67) 

s 



Consider 



dY d / n\ h n 
lit ~ di 



(7) - 7 - 



where dot means time derivative. In absence of entropy production, S — R^s is constant {R 
is the scale factor). Differentiating with respect to time we see that 

s = -3-s = -3Hs (69) 
R 

which yields 

Y=- + 3H-. (70) 

s s 

Hence we can rewrite Eq. as 

Y - -s(aeff«) (r2 - Y,l) . (71) 

The right-hand side depends only on temperature, and it is therefore convenient to use 
temperature T instead of time t as independent variable. Defining x = mi/T we have 

dY mi 1 ds , , /, ,9 ,,9N 
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where we have used 

1 1 rfs 1 ds 



(73) 



T sdT 3Hs dT 

which follows from Eq. (p9|). With the Friedmann equation in a radiation dominated universe 



- (74) 

where G is the gravitational constant, and the usual parameterization of the energy and 
entropy densities in terms of the effective degrees of freedom g^s and h^s, 

P = 9MT)^T\ s = KsiT)^T^ (75) 
we can cast Eq. (^2|) into the form [|lOj 

where Yeq can be written as 



using Eqs. (|4lJ), (|63) and ([75|). 

The parameter 17^^ is defined as 



For (Jeff, /loff and we use the values in Ref. po[| with a QCD phase-transition temperature 
Tqcd = 150 MeV. Our results are insensitive to the value of Tqcd, because due to a lower 
limit on the ncutralino mass the ncutralino freeze-out temperature is always much larger 
than Tqcd- 

To obtain the relic density we integrate Eq. ( fldj) from a; = to xq = m^/To where To is the 
photon temperature of the Universe today. The relic density today in units of the critical 
density is then given by 

= pV Peril = m^SoYo/ Pcrit (79) 

where pcrit = 3H^/8ttG is the critical density, sq is the entropy density today and Yq is the 
result of the integration of Eq. ([76[). With a background radiation temperature of Tq = 2.726 
K we finally obtain 

n^h^ = 2.755 X 10* ^^^0- (80) 
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Table 1: All Feynman diagrams for which we calculate the annihilation cross section. s{x), 
t{x) and u{x) denote a tree-level Feynman diagram in which particle x is exchanged in the 
S-, t- and u-channel respectively. Indices k run from 1 to 4, and indices c, d, e from 1 to 
2. u, u, d, d, i^, I', ^, £, / and / are generic notations for up-type quarks, up-type squarks, 
down-type quarks, down-type squarks, neutrinos, sneutrinos, leptons, sleptons, fermions and 
sfermions. A sum of diagrams over (s)fermion generation indices and over the neutralino and 
chargino indices k and e is understood (no sum over indices c, d). 



14 



4 Annihilation cross sections 

We have calculated all two-body final state cross sections at tree level for neutralino-ncutral- 
ino, neutralino-chargino and chargino-chargino annihilation. A complete list is given in Ta- 



Since we have so many different diagrams contributing, we have to use some method where 
the diagrams can be calculated efficiently. To achieve this, we classify diagrams according to 
their topology (s-, t- or u-channel) and to the spin of the particles involved. We then compute 
the helicity amplitudes for each type of diagram analytically with Reduce [[17| using general 
expressions for the vertex couplings. Further details will be found in Ref. ]18| . 

The strength of the helicity amplitude method is that the analytical calculation of a given 
type of diagram has to be performed only once and the sum of the contributing diagrams for 
each set of initial and final states can be done numerically afterwards. 

5 Numerical methods 

In this section we describe the numerical methods we use to evaluate the effective invariant 
rate and its thermal average, and to integrate the density evolution equation. 

We obtain the effective invariant rate numerically as follows. We generate Fortran routines 
for the helicity amplitudes of all types of diagrams automatically with Reduce, as explained 
in the previous section. We sum the Feynman diagrams numerically for each annihilation 
channel ij — > kl. We then sum the squares of the helicity amplitudes so obtained, and sum 
the contributions of all annihilation channels. Explicitly, we compute 



where 9 is the angle between particles k and i. (We set gi — 2 as appropriate for a neutralino.) 

We integrate over cos 9 numerically by means of adaptive gaussian integration. In rare cases, 
we find resonances in the t- or w-channels. For the process ij — > kl, this can occur when 
rrii < TUk and rrij > mi or rrii < mi and mj > m^'- at certain values of cos0, the momentum 
transfer is time-like and matches the mass of the exchanged particle. We have regulated the 
divergence by assigning a small width of a few GeV to the ncutralinos and charginos. Our 
results are not sensitive to the choice of this width. 

The calculation of the effective invariant rate W^s is the most time-consuming part. For- 
tunately, thanks to the remarkable feature of Eq. (|65|), Wcfi{Pcs) does not depend on the 
temperature T, and it can be tabulated once for each model. We have to make sure that the 
maximum p^s in the table is large enough to include all important resonances, thresholds 
and coannihilation thresholds. In the thermal average, the effective invariant rate is weighted 
by Kip^g (see Eq. (|65|)). The fast exponential decay of Ki at high pcff Boltzmann suppresses 
resonances and thresholds, as we have already seen in the example in Fig. ||. With a typical 
freeze-out temperature T — m^/20, contributions to the thermal average from values of Peff 
beyond ~ 1.5m^ are negligible, even in the most extreme case we met in which the effective 
invariant rate at high pcs was 10^° times higher than at peff = 0. For coannihilations, this 
value of PeS corresponds to a mass of the coannihilating particle of ~ 1.8m^. To be on the 
safe side all over parameter space, we include coannihilations whenever the mass of the coan- 
nihilating particle is less than 2.1m^, even if typically coannihilations are important only for 
masses less than lAm^. For extra safety, we tabulate Wcs from poff = up to pcff = 20™^,-, 



ble|[ 
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(81) 
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Table 2: The ranges of parameter values in our scans of supersymmetric models. For fj, 
and M2 the scans are uniform in the logarithms of the parameters and for the rest they 
are uniform in the parameters themselves. The number of models refers to the number of 
generated models satisfying experimental constraints. 

more densely in the important low pcft region than elsewhere. We further add several points 
around resonances and thresholds, both explicitly and in an adaptive manner. 

To perform the thermal average in Eq. (|65|), we integrate over pcS by means of adaptive 
gaussian integration, using a spline to interpolate in the {pcS, Wcs) table. To avoid numerical 
problems in the integration routine or in the spline routine, we split the integration interval at 
each sharp threshold. We also explicitly check for each MSSM model that the spline routine 
behaves well at thresholds and resonances. 

We finally integrate the density evolution equation ( [76| ) numerically from x — 2, where the 
density still tracks the equilibrium density, to xq = m^^/Tp. We use an implicit trapezoidal 
method with adaptive stepsize. The method is implicit because of the stiffness of the evolution 
equation. The relic density at present is then evaluated with Eq. (|80|). 

A more detailed description of the numerical methods will be found in a future publication 

6 Selection of models 

In Section ^ we made some simplifying assumptions to reduce the number of parameters 
in the MSSM to the seven parameters /i, M2, tan/3, itia, mo, At, and At- It is however a 
non-trivial task to scan even this reduced parameter space to a high degree of completeness. 
With the goal to explore a significant fraction of parameter space, we perform many different 
scans, some general and some specialized to interesting parts of parameter space. The ranges 
of parameter values in our scans are given in Table |^. 

We perform a 'normal' scan where we let the seven free parameters above vary at random 
within wide ranges, a 'generous' scan with even more generous bounds on the parameters, a 
'light Higgs' scan where we restrict to low pseudoscalar Higgs masses, and two 'high mass' 
scans where we explore heavy neutralinos. In addition we perform two other special scans: one 
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Figure 3: Neutralino relic density including neutralino and chargino coannihilations versus 
(a) neutralino mass and (b) neutralino composition Zg/(1 — Zg). The horizontal lines 
bound the cosmologically interesting region 0.025 < fi^/i^ < 1. 



to finely sample the cosmologically interesting light higgsino region, the other to study heavy 
mixed and gaugino-like neutralinos for which we found that coannihilations are important. 

Remember, though, that the look of our figures might change if different scans were used. 
One should especially pay no attention to the density of points in different regions: it is just 
an artifact of our scanning. 

We keep only models that satisfy the experimental constraints on squark, slepton, gluino, 
chargino, neutralino, and Higgs boson masses, on the width and on the 6 — > 57 branching 
ratio [2^, |T) . The last row in Table ^ gives the number of models which pass all exper- 
imental constraints. We include the most recent LEP2 constraints of which the most 
important one is 



> 85 GeV. (82) 

This bound effectively excludes most of the higgsinos lighter than the W studied in Refs. |^. 
LEP2 also puts a new constraint on the lightest Higgs boson mass. 



ruHO > 62.5 GeV, (83) 

valid for all a and /3. This constraint could be made more stringent if allowed to depend on 
sin^(/3— a), but we do not include this more refined version because in this study we are not 
very sensitive to this constraint. 



7 Results 



We now present the results of our relic density calculations for all the models in Table g. 
This is the first detailed evaluation of the neutralino relic density including neutralino and 
chargino coannihilations for general neutralino masses and compositions. So we focus on the 
effect of coannihilations. 
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Figure 4: Ratio of the neutralino relic densities with and witliout neutrahno and chargino 
coannihilations versus (a) neutrahno mass and (b) neutralino composition Zg/[1 — Zg). 



Fundamentally, we are interested in how the inclusion of coannihilations modifies the cosmo- 
logically interesting region and the cosmological bounds on the neutralino mass. We define 
the cosmologically interesting region as 0.025 < fi^ft,^ < 1. In this range of ^x^^ ^eu- 
tralino can constitute most of the dark matter in galaxies and the age of the Universe is long 
enough to be compatible with observations. The lower bound of 0.025 is somewhat arbitrary, 
and even if ^^.^"^ would be less than 0.025 the neutralinos would still be relic particles, but 
only a minor fraction of the dark matter in the Universe. 

We start with a short general discussion and then present more details in the following 
subsections. 

Fig. ^ shows the neutralino relic density Vl^h? with coannihilations included versus the neu- 
tralino mass TTix and the neutralino composition Zg/ {\ — Zg). The lower edge on neutralino 
masses comes essentially from the LEP bound on the chargino mass, Eq. (^). The few 
scattered points at the smallest masses have low tan/3. The bands and holes in the point dis- 
tributions, and the lower edge in — Zg), are mere artifacts of our sampling in parameter 
space. 

The neutralino is a good dark matter candidate in the region limited by the two horizontal 
lines (the cosmologically interesting region). There are clearly models with cosmologically 
interesting relic densities for a wide range of neutralino masses (up to 7 TeV) and compositions 
(up to 10"'' in higgsino fraction Zh = \ — Zg). A plot of the cosmologically interesting region 



in the neutralino mass-composition plane is in subsection 7.5 



The effect of neutralino and chargino coannihilations on the value of the relic density is 
summarized in Fig. ^ where we plot the ratio of the neutralino relic densities with and without 
coannihilations versus the neutralino mass and the neutralino composition Zg/{1 — Zg). 
In many models, coannihilations reduce the relic density by more than a factor of ten, and in 
some others they increase it by a small factor. Coannihilations increase the relic density if the 
effective annihilation cross section (acsv) < (gufn). Recalling that {cr^ffv) is the average of 
the coannihilation cross sections (see Eq. (|40|)), this occurs when most of the coannihilation 
cross sections are smaller than (cruvn) and the mass differences are small. 

Table lists some representative models where coannihilations are important, one (or two) 
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light heavy IfJ-l ^ l^^il ImI S> |Afi! gaugino 

higgsino higgsino bino 



Example No. 1 2 3 4 5 6 

At [GeV] 777 10241 358?7 4T4?f -7776.7 -1711.1 

Mz [GeV] -441.4 3894.1 -691.1 -1154.6 133.5 396.6 

tan/3 1.31 40.0 2.00 7.30 37.0 22.8 

ruA [GeV] 656.8 737.2 577.7 828.9 2039.5 435.1 

mo [GeV] 610.8 1348.3 1080.9 2237.9 4698.0 2771.6 

Ai/mo -1.77 -1.53 -1.03 -1.26 0.46 1.97 

At /mo 2.75 -2.01 -2.77 -0.80 0.11 0.52 

m^o [GeV] 76.3 1020.8 340.2 407.8 67.2 199.5 

Zg ' 0.00160 0.00155 0.651 0.0262 0.999968 0.99933 

m^o [GeV] 96.3 1026.4 364.5 418.2 133.5 396.0 

[GeV] 89.2 1023.7 362.2 414.1 133.5 396.0 

Q^h^ (no coann.) 0.178 0.130 0.158 0.00522 1.33 x 10'' 0.418 

D,^h^ 0.0299 0.0388 0.0528 0.00905 1.15 x 10* 0.418 



Table 3: Some representative models for which coannihilations are important (examples 1- 
5) and one model (example 6) for which they are not. We give the seven model parameters, 
the masses of the lightest neutralinos and of the lightest chargino, the gaugino fraction of the 
lightest neutralino and the relic densities with and without coannihilations. 

for each case described in the following subsections, plus one model where coannihilations are 
negligible. Example 1 contains a light higgsino-like neutralino, example 2 a heavy higgsino- 
like neutralino. Examples 3 and 4 have ~ \Mi \ , and example 5 has a very pure gaugino-like 
neutralino. Example 6 is a model with a gaugino-like neutralino for which coannihilations 
are not important. 

We have looked for a simple general criterion for when coannihilations should be included, 
one that goes beyond the trivial statement of an almost degeneracy in mass between the 
lightest neutralino and other supersymmetric particles. We have only found few rules of 
thumb, each with important exceptions. We give here the best two. 

The first rule of thumb is that when coannihilations are important, |///Afi| < 2. But excep- 
tions are found, as can be seen in Fig. ^ where we show the reduction in relic density due 
to the inclusion of coannihilations as a function of \fj,/Mi\. Notice that when \^/Mi\ <C 1, 
the neutralino is higgsino-like; when \iji/Mi\ ^ 1, the neutralino is gaugino-like; and when 
|/i/Mi| ~ 1, the neutralino can be higgsino-like, gaugino-like or mixed. 

The second rule of thumb is that coannihilations are important when Zg < 0.23 for < 200 
GeV and when Zg/{1 - Zg) < GeYf for > 200 GeV. There are exceptions 

to this rule, as can be seen in Fig. ^ where the ratio of relic densities with and without 
coannihilations is plotted versus the neutralino mass, the left panel for points satisfying the 
present criterion, the right panel for those not satisfying it. 

In the following subsections, we present the cases where we found that coannihilations are 
important and explain why. We first discuss the already known case of light higgsino-like neu- 
tralinos, continue with heavier higgsino-like neutralinos, the case \fj,\ ~ \Mi \ and finally very 
pure gaugino-like neutralinos. We then end this section by a discussion of the cosmologically 
interesting region. 
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Figure 5: Ratio of the relic densities with and without coannihilations versus |/i/Mi|. 
Coannihilations are important when < 2. 



7.1 Light higgsino-like neutralinos 

We first discuss light higgsino-likc neutralinos, < m^, Zg < 0.01, since coannihilation 
processes for these have been investigated earlier by other authors ||^, |^, |^ . 

Mizuta and Yamaguchi stressed the great importance of including coannihilations for 
higgsinos lighter than the W boson. For these light higgsinos, neutralino-neutralino annihi- 
lation into fermions is strongly suppressed whereas chargino-neutralino and chargino-chargino 
annihilations into fermions are not. Since the masses of the lightest neutralino and the light- 
est chargino are of the same order, the relic density is greatly reduced when coannihilations 
are included. Mizuta and Yamaguchi claim that because of this reduction light higgsinos are 
cosmologically of no interest. 

Drees and Nojiri |^ included coannihilations between the lightest and next-to-lightest neu- 
tralino, but overlooked those between the lightest neutralino and chargino, which are always 
more important. In spite of this, they concluded that the relic density of a higgsino-like 
neutralino will always be uninterestingly small unless > 500 GeV or so. 

Drees at al. |^ then re-investigated the relic density of light higgsino-like neutralinos. They 
found that light higgsinos could have relic densities as high as 0.2, and so be cosmologically 
interesting, provided one-loop corrections to the neutralino masses are included. 

We agree with these papers qualitatively, but we reach different conclusions. We show our 
results in Fig. |^, where we plot the relic density of higgsino-like neutralinos versus their 
mass with coannihilations included, as well as the ratio between the relic densities with 
and without coannihilations. The Mizuta and Yamaguchi reduction can be seen in Fig. Mb 
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Figure 6: Ratio of the relic densities with and without coannihilations versus neutralino 
mass m-)^. Coannihilations are generally not important when Zg > f{m-),), where /(nix) is 
the 'second rule of thumb' given in the text. 



below 100 GeV, but due to the recent LEP2 bound on the chargino mass the effect is not 
as dramatic as it was for them. If for the sake of comparison we relax the LEP2 bound, the 
reduction continues down to 10^^ at lower higgsino masses and we confirm qualitatively the 
Mizuta and Yamaguchi conclusion — coannihilations are very important for light higgsinos 
— but we differ from them quantitatively since we find models in which light higgsinos have 
a cosmologically interesting relic density. For the specific light higgsino models in Drees et 
al. 1^ we agree on the relic density to within 20-30%. We find however other light higgsino- 
like models with higher fl^h'^ 0.3, even without including the loop corrections to the 
neutralino masses. 

So there is a window of light higgsino models, ^ 75 GeV, that are cosmologically in- 
teresting. All these models have tan/3 < 1.6 and those with the highest relic densities have 
tan/3 ~ 1.2. These models escape the LEP2 bound on the chargino mass, ~ 85 GeV, 
because for tan f3 <2 the mass of the lightest neutralino can be lower than the mass of the 
lightest chargino by tens of GeV. By the same token, coannihilation processes are not so 
important and the relic density in these models remains cosmologically interesting. Most of 
these models will be probed in the near future when LEP2 runs at higher energies, but some 
have too large a chargino mass (mj > 95 GeV) and too large an TJj boson mass {rriffo > 90 
GeV) to be tested at LEP2. Thus ~ 75 GeV higgsinos with tan/3 < 2 may remain good dark 
matter candidates even after LEP2. 
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Figure 7: For higgsino-like neutralinos {Zg < 0.01) we show (a) the relic density with coan- 
nihilations included and (b) the ratio of the rehc densities with and without coannihilations 
versus the neutralino mass. The horizontal hues in (a) limit the cosmologicahy interesting 
region 0.025 < O^/i^ < 1. 



7.2 Heavy higgsino-like neutralinos 

Coannihilations for higgsino-like neutralinos heavier than the W boson have been mentioned 
by Drees and Nojiri , who argued that they should not change the relic density by much, and 
by McDonald, Ohve and Srednicki |6j, who warn that they might change it by an estimated 
factor of 2. We are the first to present a quantitative evaluation in this case. We typically 
find a decrease by factors of 2-5, and in some models even by a factor of 10 (see the right 
hand part of Fig. ^). 

For rrix > mw , the lightest and next-to- lightest neutralinos and the lightest chargino are close 
in mass, and they annihilate into M^-bosons besides fermion pairs. While the annihilation 
and coannihilation cross sections into W^-pairs are comparable, the coannihilation of XiX2j 
Xixt ^^'^ X2X1 into fermion pairs is stronger than the XiXi ~^ ff annihilation. This gives 
the increase in the effective annihilation rate that we observe. 

As a result, the smallest and highest masses for which higgsino-like neutralinos heavier than 
the W boson are good dark matter candidates shift up from 300 to 450 GeV and from 3 to 
7 TeV respectively. 

Together with the result in the previous subsection, we conclude that higgsino-like neutralinos 
{Zg < 0.01) can be good dark matter candidates for masses in the ranges 60-85 GeV and 
450-7000 GeV. 



7.3 Models with ~ \Mi\ 

Coannihilations for mixed or gaugino-like neutralinos have not been included in earlier cal- 
culations. It has been believed that they are not very important in these cases. On the 
contrary, when \fi\ ~ |Mi| and > niw, there is a very pronounced mass degeneracy 
among the three lightest neutralinos and the lightest chargino. The ensuing coannihilations 
can decrease the relic density by up to two orders of magnitude or even increase it by up to 
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Figure 8: For neutralinos with 0.8 < |/i/Afi| < 1.2 we show (a) the rehc density with coan- 
nihilations included and (b) the ratio of the rehc densities with and without coannihilations 
versus the neutrahno mass. The horizontal lines in (a) limit the cosmologically interesting 
region 0.025 < O^/i^ < 1. 



a factor of 3. This is easily seen in Fig. || as the vertical strip at \^/Mi\ ~ 1. In Fig. jsjthe 
relic density including coannihilations and the ratio of the relic density with to that without 
coannihilations are shown versus the neutralino mass for models with 0.8 < |/j,/Mi| < 1.2. 

We recall that in models with ^ \Mi the lightest neutralino can be higgsino-like, mixed 
or gaugino-like. If the lightest neutralino is mixed [Zg ~ 0.5), coannihilations can increase 
the relic density, whereas if it is more higgsino-like or gaugino-like they will decrease it. This 
because the annihilation cross section for mixed neutralinos is generally higher than those 
for higgsino-like or gaugino-like neutralinos. 

The largest decrease we see for this kind of models is when \Mi\ is slightly less than and 
both are in the TeV region. In this case, the lightest neutralino is a very pure bino, and its 
annihilation cross section is very suppressed since it couples neither to the Z nor to the W 
boson. The chargino and other neutralinos close in mass have much higher annihilation cross 
sections, and thus coannihilations between them greatly reduce the relic density. This big 
reduction suffices to lower fi^/i^ to cosmologically acceptable levels if Zg < 0.96. This reduc- 
tion does not occur for masses much lower than a TeV, because the terms in the neutralino 
mass matrix proportional to the W mass prevent such pure bino states and such severe mass 
degeneracy. 

To conclude, when ~ coannihilations are very important no matter if the neutralino 

is higgsino-like, mixed or gaugino-like. The relic density can be cosmologically interesting 
for these models as long as the gaugino fraction Zg < 0.96: these neutralinos are good dark 
matter candidates. 



7.4 Gaugino-like neutralinos with ^ \Mi\ 

When \fi\ ^ \Mi\, the lightest neutralino is a very pure gaugino. According to the GUT 
relation Eq. (^, the supersymmetric particles next in mass, the next-to- lightest neutralino 
and the lightest chargino, are twice as heavy. So we expect that coannihilations between 
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Figure 9: Neutralino masses and compositions Zg/(1 — Zg) for cosmologically interesting 
models (a) with and (b) without inclusion of coannlhilations. 

them are of no importance.^ In fact, as discussed in section coannlhilations would need 
to increase the effective cross section by several orders of magnitude for these large mass 
differences. 

This actually happens in some cases, the small spread at \fi/Mi\ ~ 130 in Fig. |^. In these 
models, the lightest neutralino is a very pure bino {Zg > 0.999) and the squarks are heavy. 
Its annihilation to ferniions is suppressed by the heavy squark mass, and its annihilation to Z 
and W bosons is either kinematically forbidden or extremely suppressed because a pure bino 
does not couple to Z and W bosons. On the other side, the lightest chargino is a very pure 
wino, which annihilates to gauge bosons and fermions very efficiently. The huge increase in 
the effective cross section, compensated by the large mass difference, reduces the relic density 
by 10-20%. However, the relic density before introducing coannihilations was of the order of 
lO'^-lO'*, and this small reduction is not enough to render these special cases cosmologically 
interesting. 

7.5 Cosmologically interesting region 

We now summarize when the neutralino is a good dark matter candidate. Fig. ^ shows the 
cosmologically interesting region 0.025 < il^h^ < 1 in the neutralino mass-composition plane 
Zg/{1 — Zg) versus m^. 

The light higgsino-like region does not extend to the left and down because of the LEP2 
bound on the chargino mass. The lower edge in gaugino fraction at Zg ~ 10~^ is the border 
of our survey (how high \M2\ is allowed to be). The upper limit on Zg and the upper limit 
on the neutralino mass come from the requirement fl^h^ < 1. The hole for higgsino-like 
neutralinos with masses 85-450 GeV comes from the requirement fl^^h^ > 0.025. 

We see that coannihilations change the cosmologically interesting region in the following 
aspects: the region of light higgsino-like neutralinos is slightly reduced and the big region of 
heavier higgsinos is shifted to higher masses, the lower boundary shifting from 300 GeV to 

5ln models with non-universal gaugino masses, the lightest gaugino-likc neutralino can be almost degen- 
erate with the lightest chargino, and coannihilations can bo important, as examined e.g. in Ref. bd]. 
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450 GeV and the upper boundary from 3 TeV to 7 TeV. 

The fuzzy edge at the highest masses is due to models in which the squarks are close in 
mass to the lightest neutralino, in which case t- and u-channel squark exchange enhances 
the annihilation cross section. In these rather accidental cases, coannihilations with squarks 
are expected to be important and enhance the effective cross section even further. Thus, the 
upper bound of 7 TeV on the neutralino mass may be an underestimate. 

8 Conclusions 

We have performed a detailed evaluation of the relic density of the lightest neutralino, includ- 
ing for the first time all two-body coannihilation processes between neutralinos and charginos 
for general neutralino masses and compositions. 

We have generalized the relativistic formalism of Gondolo and Gelmini to properly treat 
(sub)threshold and resonant annihilations also in presence of coannihilations. We have found 
that coannihilations can formally be considered as thresholds in a suitably defined Lorentz- 
invariant effective annihilation rate. 

Our results confirm qualitatively the conclusion of Mizuta and Yamaguchi |7| : the inclusion 
of coannihilations when < mw is very important when the neutralino is higgsino-like. 
In contrast to their calculation we do however find a window of cosmologically interesting 
higgsino-like neutralinos where the masses are ^ 75 GeV and tan/3 < 1.6. This is 
due primarily to a milder mass degeneracy at low tan/3, and secondarily to the one-loop 
corrections to the neutralino masses pointed out in Ref. 

We also find that coannihilations arc important for heavy higgsino-like neutralinos, > 
mw, for which the relic density can decrease by typically a factor of 2-5, but sometimes even 
by a factor of 10. Higgsino-like neutralinos with > 450 GeV can have ft^h'^ > 0.025 and 
hence make up at least a major part of the dark matter in galaxies. 

When \fi\ ~ \Mi\, coannihilations will always be important: they can decrease the relic 
density by up to a factor of 100 or even increase it by up to a factor of 3. In these models, 
the neutralino is either higgsino-like, mixed or gaugino-like, and when the gaugino fraction 
Zg < 0.96, the relic density can be cosmologically interesting. 

Coannihilations between neutralinos and charginos increase the cosmological upper limit on 
the neutralino mass from 3 to 7 TeV. Coannihilations with squarks might increase it further. 

Coannihilation processes must be included for a correct evaluation of the neutralino relic 
density when » |Mi| and when |^| < 2|Mi|. In the first case, the neutrahno is a very 
pure gaugino and its relic density overdoses the Universe. In the second case, the neutralino is 
either higgsino-like, mixed or gaugino-like, and for each of these types there are many models 
where it is a good dark matter candidate. To establish this, the inclusion of coannihilations 
has been essential. 
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